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Abstract 

We present the GRj unction package which allows boundary surfaces and 
thin-shells in general relativity to be studied with a computer algebra sys- 
tem. Implementing the Darmois-Israel thin shell formalism requires a careful 
selection of definitions and algorithms to ensure that results are generated in 
a straight-forward way. We have used the package to correctly reproduce a 
wide variety of examples from the literature. We present several of these ver- 
ifications as a means of demonstrating the packages capabilities. We then use 
GRj unction to perform a new calculation - joining two Kerr solutions with 
differing masses and angular momenta along a thin shell in the slow rotation 
limit. 
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1 Introduction 



The Darmois-Israel junction/thin-shell formalism has found wide apphcation in 
general relativity and cosmology ||l|, ^. The junction of dust to Schwarzschild 
by Oppenheimer and Snyder allowed the first insights into the nature of grav- 
itational collapse to a black hole Since Israel's landmark paper |Q] the 
formalism has been applied in a number of contexts ranging from further stud- 
ies of gravitational collapse to the evolution of bubbles and domain walls in a 
cosmological setting. 

In this paper we describe the GRjunction package we have developed to 
assist relativists in the evaluation of junction conditions and the parameters 
associated with thin-shells (the package is available free of charge Q). At the 
present time the package deals only with non-null surfaces - although efforts to 
extend this to null shells are underway. GRjunction runs under the computer 
algebra system Maple ^] in conjunction with GRTensorll j^. Our goal in creat- 
ing the package was to ensure that it could easily recover all the standard shell 
results in the literature (the bulk of which assume spherical symmetry) without 
biasing the package towards spherical symmetry in any way - allowing users to 
probe the relatively unstudied area of non-spherical shells and junctions. The 
package is necessarily interactive allowing users to manipulate results and de- 
termine the conditions for junctions. 

We begin by outlining the shell formalism to establish notation and moti- 
vate choices of algorithms which we describe in the following section. We then 
demonstrate the package by repeating some standard junction and shell calcu- 
lations. Next we present some new results relating to the study of shells around 
slowly spinning black holes. To validate the package we re-executed a number 
of the standard results in the literature. A summary of these tests appears in 
the final section. 

The intention of this paper is to describe the junction package we have de- 
veloped and not to review the vast literature on junctions and thin shells. The 
references we have chosen are not always the first or simplest treatment of a 
problem and in some cases we have deliberately selected examples which differ 
from the standard treatments to test the robustness of our package. 



2 The Formalism 

In this section we review the junction formalism to establish notation. For more 



detailed discussions see e.g.M]-[10|. 



Consider two spacetimes (Lorentzian manifolds with signature (+ + H — )) 
M+ and M~ with metrics g'^p{x\) and g^jsi^"^) iii the coordinate systems x'\. 
and xZ. Within these spacetimes define two non-null 3-surfaces S"^ and S~ (in 
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and M~ respectively) with metrics (7^(^!^) and g^j{^'L) in the coordinates 
and which decompose each of the 4-manifolds into two distinct parts. 
(Greek indices range over the coordinates of the 4-manifold and Roman indices 
over the coordinates of the 3-surfaces). We label the distinct parts of M+ cre- 
ated by S+ as Ad^ and M2 and likewise for M~ . The junction/shell formalism 
constructs a new manifold Ai by joining one of the distinct parts of M'^ to one 
of the distinct parts of M~ by the identification S+ = T,^ = 11. Clearly there 
are four possibilities, i.e. M^*" U M{,M^ U Mf , M^*" U M2,M^ U ■ The 
assumed isometry between the points on the surface is often (but need not al- 
ways be) the identification S,^ = S,*!. 

What now follows holds simultaneously for and M~ and so we drop 
the lb distinction in this paragraph. The parametric equation for E is of the 
form 

/(x"(r)) = 0. (1) 

We assume that S is non-null. The unit 4-normals to S in M are given by 

1 df 

We assume Ua ^ and label S as timelike (spacelike) for A = —nan" = —1(1). 
The three basis vectors tangent to S are 



which give the induced metric on E by 



efa) = 7^ (3) 



3ij = -Q^-Q^9af3- (4) 

The extrinsic curvature (second fundamental form) is given by 

( d'^x^ dx^ \ 

We define the four tangent to E and tt"(= u^^V pu"^) the four-acceleration. 
Combining these with the above relations gives 

Kit"] = -v^v? [Kij] (7) 

and 

UaVP" = -u'u^^ (8) 

where [X] = X+ |s -X" |s and X = {X+ |s +X- |s)/2 with X± |s denoting 
the limiting values of X on E. 
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The Darmois conditions for the joining of a part of to a part of M are 

[9ij\ = (9) 

and 

[K,,] = (10) 

If both (I) and (0) are satisfied we refer to S as a boundary surface. If 
only (P) is satisfied then we refer to S thin-shell. 



Conditions (^) and (|l^) require a common coordinate system on S and this 
is easily done if one can set = . Failing this, estabhshing (P) requires a 
solution to the three dimensional metric equivalence problem. Condition (p!c|) 
as it stands is ambiguous since the orientation of the 4- vector field = | s 
has not been specified. The Israel formalism requires the normals in M to point 
from to (where A denotes the part of M~ and B denotes the part of 
we wish to use to form ^A). Clearly the sign of the normal vectors are 
crucial since e.g. n~ points away from the portion of M~ which will be used 
in forming M. Hence an understanding of which side n~ points into is key. In 
general this can be done by considering a trajectory in M~ through E~ with 
tangent n~ (and likewise for n+). The majority of the existing literature deals 
with spherical symmetry where the direction of the normal is clear, but in more 
complicated examples (see below) great care must be taken. Note that while 
there are two normal vectors in each of once we have identified S~ and S"*" 
there is a single unique normal field to S in A^. There may be circumstances 
in which one can determine the differential relation between the coordinates of 



and via 



(11) 



in an open neigbourhood of S in and this will give the direction of one 
normal vector relative to the other. (In the case of a boundary surface often 
only the sign of on one side of S needs to be determined. Since V^n" = Kl 
in the case of a boundary surface ( p!o| ) gives the useful relation [V^n"] = 0.) 

Some studies have left the sign of the normal vectors unspecified to ex- 
haustively study the taxonomies of all possible combinations of M~ and 
(usually excluding those which require shells which violate energy conditions 
e.g. |l^, |l^). The junction package allows the user to leave the sign unspecified 
but we take the view that the "typical" starting point is to explicitly choose 
the signs of the normal vectors so that a particular combination of and 
Mj^ can be studied. 



Once we have selected signs in (g) there is no ambiguity in (10) and we 
use (^ and (|lO|) in conjunction with the Einstein tensor Gap to obtain the 
identities 

Gapn'^n^] = (12) 
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and 



0. 



(13) 



(Note that (|^) and (|T^) do not guarantee (^) and (|T^).) This shows, for ex- 
ample, that for timelike S the flux through S (as measured comoving with S) 
is continuous. 



The Israel formulation of thin shells follows from the Lanczos equation 

A 



5, 



I J — g^d^ij] 9ij 



) 



(14) 



and we refer to Sij as the surface stress-energy tensor of S. The "ADM" 
constraint 



along with Einstein's equations then gives the conservation identity 



The "Hamiltonian" constraint 

G^pn'^n^ = iA{^R) + - KijK'^)/2 
gives the evolution identity 



(15) 



(16) 



(17) 



(18) 



The identities (16) and (p!q) do not give information about the dynamics of 
the shell. The evolution of the shell stems from a phenomenological interpre- 
tation and the Lanczos equation (0). 



2.1 Phenomenology 

The standard phenomenology associated with T, is introduced as follows: let 
be the coordinates intrinsic to S and consider a trajectory ^"(t) with associated 
3-tangent 

u'' = ^ (19) 
dr 

where u°'Ua = A. For a timelike surface we view the curve as the worldline of 
a (possibly hypothetical) particle in the surface with r the proper time. (If the 
surface is spacelike the phenomenology is a formal analogy). The associated 
four tangent is 

dx4_ „ dx'l , , 

^ dr ' ' 



Equivalently 



ut'-^ (21) 
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We view E as covered by a 3- vector field u"". Note that r is defined curve by 
curve on S but not in general over the entire surface. If r is defined over the 
entire surface then we label this as a "one-parameter surface". The surface 
energy density a associated with the trajectory on S is defined by 



(22) 



where Uaq"" = 0. The inclusion of qa represents an intrinsic energy fiux orthog- 
onal to the prescribed intrinsic velocity field. The surface energy density cr(^'^) 
is not in general an eigenvalue but is given by 



(23) 



In analogy to a perfect 4-fluid we suppose that Qa = and that Sab takes 
the form 

Sab = -A(a(r ) + p{e))uanb + pi^gab- (24) 



This defines the surface pressure ( - surface tension) 

p(e^) = (a(a + S:)/2. 
It now follows from (^) that 



& + (a + p)^ = -A 



Tap 



-n' 



(25) 



(26) 



where ' signifies the intrinsic 3-derivative along u"" ( & = u^'VaO') and ^ is the 
three-expansion of u"" (<I> = Vau"'). Prom (^), the Lanczos equation (|l4|) and 
the phenomenological equations (p3|) and (|25|) we have 



(27) 



[n„n"] = 8iT{p + a/2). 
From (|6|), (|l^) and the identity ([l8|) the phenomenology gives 

-A 



A{a + p)nau'^ + pKl 



(28) 



In spherical spacetimes the first integral of the evolution equation (28) is 



given by the identity |11| 



A+(l^ 
' M 



2Am 



R 



+ 



2r) 



(29) 



where R = R{t), ' = d/dr and M = AiraR? . The effective mass, m± is given 

by 

m± - \{^'^9fef'' (30) 
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3 Computer Implementation 



Specifying the junction formalism for a computer algebra system requires a 
careful choice of specific object definitions and recognition of several standard 
calculus manipulations. Joining spacetimes requires the simultaneous consider- 
ation of four metrics, so the underlying general relativity software must permit 
this. In order to make the package reproduce the existing literature on spherical 
shells special consideration must be given to the choice of object definitions for 
one-parameter shells. We first present an overview of the package and then dis- 
cuss the choices we have made in implementing the junction/thin shell package 
for the Maple version of GRTensorll. 



3.1 Overview 

The GRjunction package provides the user with a means to specify a sur- 
face, calculate intrinsic and extrinsic quantities on the surface, identify two 
such surfaces and evaluate whether a boundary surface or thin-shell results. 
The specification of a surface is done by invoking the command surf which 
then prompts the user for the necessary information. Once a surface has been 
specified objects defined on the surface (e.g. Kij) can be calculated. The iden- 
tification of two surfaces is performed via the command join which calculates 
[gij] and displays the result. If [gij] ^ the user can manipulate this expression 
in an attempt to determine restrictions on metric functions which will ensure 
[gij] = (see the first example below). The jump or mean of any quantity in 
the joined manifold can be evaluated by means of the operators Jump and Mean. 
Hence to determine [Kij] the user would refer to the object Jump [K(dn,dn)] ; 
to determine K''^ , Mean [K (up, up)] . Standard quantities and equations for the 
joined spacetimes can be calculated in a straight-forward manner. For example 
to determine Si the user refers to the object S3(dn,up). (By convention we 
list dn indices ahead of up indices for mixed two index objects). 

While we have emphasised the Darmois-Israel formalism, GRjunction can 
trivially evaluate the Lichnerowicz junction conditions |jl^ which consist of 

[9af3] = (31) 

[dgap/dx^] = (32) 

(Recall these conditions require admissible coordinates and consequently are 
not as general the Darmois-Israel conditions). The object dgap/dx'^ is re- 
ferred to as g(dn,dn,pdn) (pdn denoting the partial derivative index) and 
so (|3^ ) can be evaluated by referring to the object Jump [g(dn,dn, pdn)] . 
GRTensorll also allows users to work within the Newman- Penrose formalism 
and the GRjunction will allow users to evaluate jumps in the spin coefficients 
across a surface. (This technique is employed in e.g. [|15[). 
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3.2 The surface and related quantities 

In this section we describe how to specify a non-null surface to the junction 
package and the intrinsic (on S) and extrinsic quantities which can be calcu- 
lated with the package. In this section we restrict attention to quantities which 
are independent of phenomenology. 

All the expressions in this section should implicitly carry a it designation 
(which we omit) with the exception of the i^*, the coordinates on S. For the 
package to compare first and second fundamental forms we must have ^!|_ = .^L 
The package does not currently consider the three-metric equivalence problem. 

In general to specify a surface S in a space M with coordinates in 
sufficient detail to allow calculation of the first and second fundamental forms 
we require 

• the coordinates on S 

• the coordinate definition of S: = x"(^*). 

• the parametric definition of S: /(x") = 

• the choice of normal vector sign in (^) 

The essential idea of the Darmois/Israel formalism is to use intrinsic quan- 
tities in the description of all objects of interest on S. However many of the 
definitions of objects on the surface are in terms of quantities in M. For exam- 
ple (^ uses the normal vector which involves partial derivatives with respect 
to the coordinates of M and Christoffel symbols of M. This can frequently 
be resolved simply by substituting the coordinate definition of the surface (i.e. 
x" = x"(^*)) but this is not always desirable. Consider a metric which has an 
arbitrary function u{r,9) and a definition of T, with coordinates = {9,(j),T) 
and which has a coordinate definition which includes r = f{C) and 9 = 9. 
The Christoffel symbols used in defining Kij may contain partial derivatives 
of u with respect to both r and 9. Substitution oi 9 = 9 will merely change 
the variable, but the r substitution will result in 9ii(/(^*), 0)/9/(^*). This is 
more than notationally ugly - it constitutes an error in Maple; you cannot take 
a partial derivative with respect to a function. In these cases we make use 
of the Maple operator D (instead of the Maple procedure diff ). For example 
the representation of the r derivatives of /(r) and i'{r,9) on the surface are 
represented in Maple as 

^/(O l.=«(t) - D{f){R{t)) (33) 

^K^,^)lr.=ii(0 - D[l]{u){R{t),9) (34) 

(35) 

(the [1] indicates that the derivative acts on the first argument of the function 
I'). It is essential that we make use of this form of representing derivatives to 
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ensure that their time dependence is handled properly in subsequent calcula- 
tions. While we could use the D derivative throughout we prefer to minimize 
it's use since in the Maple output we find it preferable to have expressions 
with df{r)/dr instead of D(f ) (r). Consequently object defintions which take 
derivatives and then evaluate values on the surface make use of code which sub- 
stitutes in the coordinate relations and ensures that the substitution does not 
cause errors in derivatives. This code then minimizes the number of D deriva- 
tives which appear by converting extraneous occurances back into Maple's usual 
diff derivative. 



3.3 One-Parameter Surfaces 

A large part of the existing shell literature has dealt with timelike spherical 
3-surfaces within spherical 4-manifolds. The majority of these analyses define 
one of the to be the proper time on the surface (what wc referred to above 
as one-parameter shells). We must take some care in choosing algorithms for 
objects in these cases since some of the routine steps in a hand calculation are 
best avoided in a computer algebra approach. All of the algorithm issues arise 
in the specification of a surface in one of so we limit our discussion to the 
specification of one surface. 

First we describe two standard timelike surfaces we will make reference 
to in our discussion: a static spherically symmetric surface and the dynamic 
counterpart within a spherically symmetric 4-manifold. In the 4-manifold we 
take coordinates {r,6,(j),t) and on the 3-surface coordinates {9,4>,t). For the 
static shell we use the coordinate definitions 

r = R, 9 = §,(!) = ^,t = Ts{T) (36) 

and a surface equation r — R = 0. For the dynamic shell we use 

r = R{T),e = 9,(l) = 4>,t = Ta{T) (37) 
and a surface equation r — R(t) = 0. 

In both cases it is conventional to ensure that the coefficient of dr^ on S is 
— 1 (hence r is the proper time on the shell). For this reason we express t as a 
function of r and then make use of the constraint Uau"' = — 1 to eliminate e.g. 
dT^{T)/dT from the quantities we calculate on S. This facility is built into the 
surf command. The user is asked if the surface has a parameter which governs 
its evolution. If the user so indicates then the package evaluates Uau" and 
asks if there is a quantity to be eliminated. An attempt is then made to solve 
the tangent constraint so that this quantity can be eliminated. The tangent 
constraint is then automatically applied during the calculation of Ua and Kij. 

We must specify an algorithmic approach to calculating quantities such as 
Kij and n° which recover the results for one-parameter shells in a straight- 
forward way. Not all object definitions are equivalent because some require 
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steps which are obvious in a hand calculation but are difficult for a computer 
algebra system. One example of this is the calculation of u°. For a one param- 
eter shell with = x°'{t) we would proceed with the definition: 

-^V,n" = u^(^ + T-s.uA (38) 



- dr dxP +^/37«^ 

= + (39) 

While in a hand calculation using the definition of a total derivative to get the 
last line is an obvious step this is not so in a computer algebra system. Once 



specific components are used in (38) and partial derivatives are taken, recogniz- 
ing terms which can collected into total derivatives with respect to r becomes a 
problem in pattern matching. To avoid such problems we require that the user 
indicate the variable to be used for total derivatives and we use the relation 
( |3g| ) instead of (^). The same issue arises in the calculation of Kij for the 
dynamic shell discussed above if we use (^) so we use the definition (^). This 
definition has the advantage that it is expressed in terms of derivatives with 
respect to intrinsic quantities (although we still need to evaluate the Christoffel 
symbols of the 4-manifold on S). 



Another computer algebra issue arises in the calculation of the normal vec- 
tor. For one-parameter surfaces it is customary to write the equation of the 
surface in terms of a function of a parameter (e.g. r — R{t) = 0). Yet the 
covariant derivative in M will require partial derivatives with respect to the 
x". The junction package requires that the surface equation be in terms of 
the (so for the dynamic shell example above we specify r — R{t) instead of 
r — R[t)). In these cases the package then applies the chain rule 

dr dR dR{t)dT dR{T)/dT 

dt^dT^ dr ~dt^ dT{T)/dT ' ^ ^' 

where we have made use of the relations r = R{t) and t = T{t) and that 
T 7^ If the surface equation is governed by a parameter other than r then the 
normal vector and Kij can be expressed in terms of intrinsic quantities only if 
the function dependence on the other ^* can be given explicitly in the surface 
equation (e.g r — r{t)cos6 = 0). If all that is known is e.g. r — R(t,6) then 
we cannot make use of ( ^0|) and and K^j will be left in terms of the partial 
derivative with respect to the coordinates of M. 

With these definitions and constraints we now have available the quantities: 
n" , , u°' , Qij and Kij in terms of intrinsic quantities. We require a few addi- 
tional objects which are calculated in a straight-forward way. A summary of 
all the objects which can be calculated on a surface by the junction package is 
provided in Tables 1 and 2. 
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Object 


Definition 


Name 


Comments 


,x"(f ) 




xf orni(up) 


coordinate definition of E 


A 




utype 




/ 


/(.x"(e)) 


surface 


parametric surface equation 




9//(?x" (normalized) 


n(dn) 


the normal to E 






u(up) 


four tangent to S 








(one parameter shells only) 


T 




totalVar 


shell parameter (e.g. proper time) 








(one parameter shells only) 






udot (up) 


four acceleration of S 






(one parameter shells only) 


m 




mass 


mass (for spherical 








one parameter shells only) 



Table 1: Objects in M± 



Object 


Definition 


Name 


Comments 


e 




x(up) 


coordinates on S 


9ij 




g(dn,dn) 


first fundamental form 


Kij 


„ 1 9^3:7 p7 Q^o^ \ 
T \di^di^ "1" ^ a/3 di^ 9$^ ) 


K(dn,dn) 


second fundamental form 


Kl 


trK 






((3)i? + K'^ + KijK'^)/2 = Ga^n^'n^ 


Ksq 






HCGeqn 


Hamiltonian constraint 






HCTeqn 


Hamiltonian constraint 






PCGeqn 


Momentum constraint 






PCTeqn 


Momentum constraint 






u3(up) 


three tangent 






(one parameter surf, only) 






u3div 


three divergence 

(one parameter surf, only) 



Table 2: Objects on E 
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Object 



Definition 



Name 



Comments 



[nQ,'U°'] = —u^u^ [Ki 
naU°^ = —u^u^Kij 
divergence (T) 
divergence (G) 
mean 
mean 



a/3 Q^t II' 

a.B pici ^ 



A 



•/3 a? 



S3(dn,dn) 

nuJumpeqn 

nuMeaneqn 

divSTeqn(dn) 

divSGeqn(dn) 

SKTeqn 

SKGeqn 



intrinsic stress 

energy tensor 

identity 

identity 

identity 

identity 

identity 

identity 



Table 3: Objects in on S (General Case) 



3.4 Objects in M 

Once the user has loaded two 4-manifolds and specified two 3-surfaces the 
command join is used to identify the two surfaces. A variety of objects relating 
to the boundary surface or shell can now be calculated (Tables 3 and 4). 

The user now has four active metrics in the GRTensorll session (S^, M^). 
After using join the default metric is S"*". The operators Jump and Mean will 
by default make reference to objects in S^. In general these operators can be 
used to take the jump or mean of objects from the default metric to any user 
specified metric. For example if the user had changed the default metric to M"*" 
(say Schw) and wished to evaluate ( |3l| ) with M~ as e.g. Schwinterior this 
would be done by using the metric name of M~ as a second parameter to Jump 
i.e. Jump [g (dn , dn) , Schwinterior] . 

In practice users may wish to structure calculations and determine results 
in a variety of ways and this can require minor differences in the definitions 
of objects. Consider a user who wishes to study the dynamics of a Minkowski 
- shell - Schwarzschild scenario. After specifying the two surfaces and joining 
them they might elect to examine the density and pressures of the shell as given 
by (p^) and ( p5[) which is done by referring to the objects sigma and P. When 
considering the dynamics they may wish to have a and P in explicit form or 
they might opt to have them appear simply as cr{£,°') and -P(^") since this al- 
lows them to set -P(C") = to study the dust case. Consequently there are two 
versions of the history and conservation equations in Table 4 (e.g. HGeqn vs 
HlGeqn). 

Further variations are provided to allow users to specify a stress-energy of 
the 4-manifold for use in the conservation and history equations. A user might 
opt to specify e.g. Tg = p{r) instead of using Gg/Svr which might be a less 
transparent combination of functions of the 4-metric. This requires the defini- 
tion of separate objects for each case (e.g. HGeqn vs HTeqn). 
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Objecl 



Dcfiuili 



Name- 



Comment f 



a 

P 

Pi 



conservation 
law (T) 
conservation 
law (G) 
conservation 
law 1 (T) 
conservation 
law 1 (G) 
history (T) 
history (G) 
history (T) 
history (G) 

evolution 
integral 

evolution 
integral 



{a + St)/2 

Ktt"] = 87r(P + c7/2) 

Kii"] = 8^(Pi + ai/2) 

a + {a + P)$ = -A [Tc^^it^n^ 



0=1 + (o-i + Pi)$ = -A 



0=1 + ((71 + Pi)$ = -A 



/Stt 



A(o- + P)no,u°' + PKl 



A(o- + + PKl 



-A 
-A 



A(oi + Pi)n„M° + PiK| 



Gapn'^n'^ /Stt 
-A \Taf^n°'n'^ 
-A G^^n^n/^l /Svr 



A((Ti + Pi)n„n« + PiK| 



2Arn 
R 



Sigma 
P 

sigmal 
PI 

nuPeqn 

nuPleqn 

CTeqn 

CGeqn 

ClTeqn 

ClGeqn 

HTeqn 

HGeqn 

HlTeqn 

HlGeqn 

evlntl 

evint 



surface 

energy density 
isotropic 
surface pressure 
arbitrary energy 
density on S 
arbitrary isotropic 
pressure on S 



first integral of evolution 

equation (spherical 

one parameter shells only) 

first integral of evolution 

equation (spherical 

one parameter shells only) 



Table 4: Objects in on E (one parameter surfaces) 
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The implementation of (^) is complicated by the fact that the definition 
requires coordinates 9 and (p a-iid makes reference to the shell func- 

tion R{t). Users must adhere to these conventions (by using coordinate names 
theta, phi and r and specifying r = R{t)) if they wish to make use of this 
relation for the first integral of evolution equation. The junction package checks 
for compliance with these conventions before it will evaluate (^9|) to ensure that 
spurious results will not be generated. 



4 Examples 

4.1 Joining Schwarzschild to Uniform Dust 

We begin the demonstration of the GRj unction package by deriving the junc- 
tion of the Schwarzschild metric to an FRW spacetime (Many standard 
treatments merely verify that given certain values for m and R{t) in Schwarzschild 
a boundary surface exists, see e.g. 0]). The example also illustrates the interac- 
tive process which is required to determine junction conditions. The pedagogic 
comments make the example appear somewhat lengthy, however the example 
is computationally trivial (less than 3 CPU seconds on a SUN SPARC 5). 

We seek to join the Schwarzschild metric in the form 
dr"^ 

ds^ = + r'^dOl + sin 02^#i - (1 - 2m/r)dt2^ (41) 

1 — 2m/r 

to a portion of the closed FRW metric written as 

ds2 ^ a{tf{dx^ + smx^idOl + sin ^i^dc/.?)) - dt^ (42) 

We take {6,(f>,T) as coordinates on S. The definition of the surface in (|4l| ) is 

r = R{T),e2 = 9,cl)2 = (l),t2=T{T) (43) 

with r — R{t) = as the parametric form for the surface. For the FRW metric 
the surface is specified as 

X = X,0i = 0,01 = 0,t = T (44) 
with parametric form x ~ ^ = ^■ 

The following session output (from this point to the end of this subsection) 
was taken directly from the Latex output provided by Maple. 
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We begin by loading the GRTensorll library into Maple. 

> readlib(grii) : 

Next we load the junction package into GRTensorll. 

> grlib (junction) ; 
The GRJunction Package 

Last modified September 28, 1995 

Developed by Peter Musgrave and Kayll Lake, (c) 1995 

The first step is to load the Schwarzschild metric via qload. 

> qload (Schw) ; 

Default spacetime = Schw 
For the Schw spacetime : 

Coordinates 

= rjx"^ = 92,X^ = (j)2,X^ = t2 
Line element 

ds^ = ^^+r^ d92' +r^smi92f d <P2' + ( -1 + 2-) dis" 
1-2— V r J 

r 

Now we use the surf command to specify a surface, surf prompts the user 
for information which defines the surface and its normal vector. 

> surf (Schw, ssurf) ; 

First we arc asked for the coordinates on the surface as a list. 
Please enter the coordinates of the surface as a list 
e.g. [theta, phi, tau] ; 
Enter a list > 

> [theta, phi , tau] ; 

Next we are asked to specify which of the is the parameter for a one-parameter 
shell (if any). 

For a one-parameter shell enter the parameter (0 for none) > 

> tau; 

Now we provide the coordinate definition of S. (Note that in Maple [ ]'s do 
double duty; they denote lists and are used to indicate subscripts) 

Please enter the coordinate definition of the surface 
(the xVa.} = x(xi-["b}) ) as a LIST, 
e.g. [ r=R(tau) , theta=theta, phi=phi, t=T(tau)] 
> 

> [r=R(tau), theta[2]=theta, phi[2]=phi, t[2]=T(tau) ]; 

CPU Time = .050 

The character of the normal vector is now entered. 
Please indicate the nature of the surface normal vector 
(-1 = timelike, 1= spacelike) Enter +1,0 or -1 > 
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> 1; 

Since we have provided a shell parameter we are given the option of employing 
the relation -Uqu" = — 1 to eliminate the derivative of one of the quantities used 
in the coordinate definition of the surface. We choose to eliminate 9T(r)/9r 
(in Maple parlance dif f (T(tau) ,tau)). This choice will produce dr^ = — 1 on 
S. 

Use +/- 1 = u{"a} u-[a} as a constraint ? Enter 1 if yes > 

> 1; 

Created definition for u(dn) 

CPU Time = .100 
The constraint equation is : 

^_ -r^ (|:R(r))' + (|:T(r))'r^-4(|:T(r))'rm + 4(|:T(r))'m^ 

{—r + 2m)r 

Enter the term you wish to use the constraint to eliminate 
Term > 

> diff (T(tau) ,tau) ; 

In solving for dT^r) /dr Maple determines that there are two choices (differing 

in sign). We are asked to choose one of them. 
Solve returned multiple solutions . They are : 



r-2m + r (|:R(r))^ 



+ 4 



m . m? 



r + 2m —r + 2m {—r + 2m)r 



r-2m + r (|:R(t))^ 



+ 4 4 



—r + 2m —r + 2m {—r + 2m)r 



Please select a solution. 
Enter choice > 

> 1; 

Now wc enter the parametric definition of the surface. 
Default metric is now Schw 
Please enter the equation for the surface. 
The surface will be defined by setting the 
equation you enter to zero. 



surface 
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> 



> r-R(t[2]); 

CPU Time = .083 

Finally, we are asked explicitly for the sign of the normal vector. 
The definition of the normal vector is +/- grad( surf ace) 
please enter +1 or -1 to indicate the CHOICE of sign 
Enter +1,-1 > 

> 1; 

CPU Time = .184 

Default metric is now ssurf 

CPU Time = .133 
For the Schw spacetime : 
The Equation of the surface 

surface = r — R (t2) 
For the Schw spacetime : 
Coordinate transforms onto the surface 
xform = R( T ) 
xform = 6 
xform '^'^ = (f) 
xform *^ = T( T ) 
For the ssurf spacetime : 
Line element 

ds'' = R{Tf de^ +R{Tfsm{ef dct>^ - dT^ 

The intrinsic metric and normal vector have been calculated. 
You may wish to simplify them further before saving the surface 
or calculating K(dn,dn) 

This completes the specification of the surface in the Schwarzschild exterior. 
Next we load the FRW metric (which we will take as the interior). 

> qload(f rw) ; 

Default spacetime = frw 
For the frw spacetime : 
Coordinates 
x^ = X,x^ = 9i,x^ = (j)i,x'^ = t 
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Line element 



ds^ = aitf dx^ +a{tfsm{x? dOi^ + &{tf s\n{xf s\ii{eif dcpi^ - dt^ 
Before defining the surface in the FRW spacetime we first calculate for 
later use. 

> grcalc(G(dn,up)) ; 

Created definition for G(dn,up) 

CPU Time = .634 

We specify the surface in M~ in parametric from as x — X = 0. The process 
is identical to that followed for the Schwarzschild case above (except we do not 
use UaU"^ = — 1). 

> surf(frw, f surf ) ; 

Please enter the coordinates of the surface as a list 
e.g. [theta, phi, tau] ; 
Enter a list > 

> [theta, phi , tau] ; 

For a one-parameter shell enter the parameter (0 for none) > 

> tau; 

Please enter the coordinate definition of the surface 

(the x{"a} = x(xi{"b}) ) as a LIST, 
e.g. [ r=R(tau) , theta=theta, phi=phi, t=T(tau)] 
> 

> [ chi = X, theta [l]=theta, phi[l]=phi, t=tau] ; 

CPU Time = .050 

Please indicate the nature of the surface normal vector 
(-1 = timelike, 1= spacelike) Enter +1,0 or -1 > 

> 1; 

Use +/- 1 = u{"a]- u{a} as a constraint ? Enter 1 if yes > 

> 0; 

Default metric is now frw 
Please enter the equation for the surface. 
The surface will be defined by setting the 
equation you enter to zero. 



surface 

> 

> chi-X; 

CPU Time = .016 
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The definition of the normal vector is +/- grad(surf ace) 
please enter +1 or -1 to indicate the CHOICE of sign 
Enter +1,-1 > 

> 1; 

CPU Time = .050 

Default metric is now fsurf 

CPU Time = .066 
For the frw spacetime : 
The Equation of the surface 
surface = x — X 
For the frw spacetime : 
Coordinate transforms onto the surface 



xform ^ 


= X 


xjorm 


= e 


xform '^^ 


= 4> 


xform * 


= T 



For the fsurf spacetime : 
Line element 

ds^ = a.{Tfsm{X f de^ +&{T fsm{X fsm{e f d^^ - dr^ 

The intrinsic metric and normal vector have been calculated. 
You may wish to simplify them further before saving the surface 
or calculating K(dn,dn) 

Now we can identify the two surfaces we have specified by means of the com- 
mand join. By convention the first surface name in the join command is 
taken as S+ for the purposes of evaluating e.g. [gij]- 

> join(ssurf , fsurf ) ; 
ssurf and fsurf are now joined. 
The default metric name is ssurf . 
The exterior metric is : ssurf 
The interior metric is: fsurf 

CPU Time = .034 

For the ssurf spacetime : 
Jump from defaultMetric — Mint 
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Jump [g{dn, dn),fsurf]ge = R( r )^ — a( r )^ sm( X )^ 

Jump [g{dn, dn),fsurf]^ ^ = K{t)'^ sin( )^ — a( r )^ sin( X )^ sm( 9 )^ 
To obtain [gij] = we require a particular value for R{t). We use the GRTensor 
command gr component to extract [gn] (which we assign to juinp_gll). We then 
set juiiip_gll equal to zero and solve for R{t). 

> jump_gll := grcomponent(Jump[g(dn,dn)] , [1,1]); 

jump.gll := R(t)2 - a(r)^sin(X)2 

> sol := [solve (jump_gll = 0,R(tau))]; 

sol := [a(T)sin( X ),— a(r)sin(X)] 

We now substitute the positive solution back in to verify [gij] = 0. This is 
accomplished by using the routine grmap to map the Maple substitution com- 
mand subs over the components of Jump [g(dn,dn)] . (The 'x' is a placeholder 
indicating which of the arguments to subs is to be filled in by the component 
value.) 

> grmap ( Jump [g(dn,dn)] , subs, R(tau) = sol[l], 'x'); 
Applying routine subs to Jump [g(dn,dn) ,f surf ] 

(Here we use a GRTensor short-cut. The _ refers to the last mentioned object. 
In this case Jump[g(dn,dn)] allowing us to save some typing). 

> grdisplay (_) ; 

For the ssurf spacetime : 

Jump from defauUMetric — Mint 

Jump[g{dn, dn),fsurf] = All components are zero 
To establish a Schwarzschild-Dust boundary surface we next need to establish 
that [Kij] = 0. We begin by calculating Jump [K(dn,up)] . We prefer the mixed 
form Kf so that we can make later use of and the phenomenology of the 
FRW space. 

> grcalc(Jump[K(dn,up)] ) ; 
Created definition for K(dn,up) 

CPU Time =1.183 

> grdisplay (_) ; 

For the ssurf spacetime : 
Jump from defauUMetric — Mint 

Jump [K{dn, up),fsurf]g ^ = 

R(r)-2m + R(T) (|:R(r))' ^ ^ ^ ^ .... 
-\ , ^ ^a(T)sin(X) + cos(X)R(T) 

_ _i ) 

R(r)a(T)sin(X) 
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Jump {K{dn, up),fsurf]^'^ = 

R(r) (^R(t)Y 

■a(T)sin(X) + cos(X)R(t; 



R(T)-2m + R(r) (|rR(T 
^ R(7) 



Jump [K(dn, up),fsurf]^^ 



R(r)a(T)sin(X) 

R(r)2 (|^R(r))+m 



R(r)-2m + R(r) (|:R(r 
^ R(7) 



R(r) 
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Setting [iiT:^] = gives a value for the Schwarzschild mass m in terms of R{t) 
and we substitute in this value. This leaves only the angular jump in K to 
contend with. 

> grmap(_, subs, m = -R(tau) '~2*dif f (R(tau) ,tau,tau) , 'x'); 
Applying routine subs to Jump [K(dn, up) ,f surf] 

Once again we make use of the value of R{t) which was required for [gij] = 0. 

> grmap(_,subs, R(tau) = sol[l], 'x'); 
Applying routine subs to Jump [K(dn, up) ,f surf] 

> gralter(_,f actor) ; 

Component Alteration of a grtensor object: 

Applying routine factor to object Jump [K(dn, up) ,f surf ] 

CPU Time = .150 

> grdisplay (_) ; 

For the ssurf spacetime : 
Jump from defaultMetric — Mint 

Jump [K{dn, up),fsurf]Q^ = 

-^l + 2a(r)sin(X)2 a( r )) + sin( X )2 (|, a(r ))' + cos( X ) 

a( r ) sin( X ) 

Jump [K(dn, up),fsurf]^ = 

-^l + 2a(r)sin(X)2 a( r )) + sin( X )2 (|, a(r ))' + cos( X ) 

a( r ) sin( X ) 

At this point we need to make reference to the Einstein tensor for the FRW 
spacetime. To this point we have not imposed the restriction that the interior 
be dust. We do this by setting 0^ = and solving for d"^ a{t) / dt"^ . This is then 
used in Jump [K(dn,up)] (with f = r on the surface as required). Recall that 
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by default object names refer to the ssurf spacetime (i.e S^). The use of a 
metric name in square brackets after the tensor name below indicates which 
metric the object is to be taken from. 
> grdisplay (G [f rw] (dn,up) ) ; 

For the frw spacetime : 

G{dn, up) 

sinix? (|a(0)' + 2a(t)sin(x)2 a( t )) + 1 - cos( x 



G ^ — — 



a( t )2 sin( X 



^ l + +2a(t) (^a(O) 



G 



a(t)2 

l+(|a(0)' + 2a(t) (^a(t)) 



PI _ 



a(02 

, 2 



^ ^_ 2sin(x)' + 3sin(x)' (|a(t)) +l-cos(x)' 
* a(0^sin(x)2 
To make use of the condition p = we will extract the G^j component (using 
grcomponent) and set it equal to zero and then isolate the d'^a{t)/dt? term. 
We will then substitute this into Kl . (We need the Maple isolate library 
which we we now load) 

> readlib(isolate) : 

The d'^a{t)/dt'^ term is now isolated and we change t to r (via the subs 
command) . 

> da := subs( t=tau, isolateC grcomponent (G [f rw] (dn, up) , [2,2] ) = 0, 

> diff (a(t) ,t$2) )); 



^2 ^ 1 -l-(^a(r 
^« ■= ^a(r) 
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dr"^ ^ ' 2 a(r) 
The resulting equation da is substituted into Kl and some simplification is 
then performed. 

> grmapC Jump [K(dn,up)] , subs, da, 'x'); 
Applying routine subs to Jump [K(dn, up) ,f surf ] 

> gralter (_ , expand , trig) ; 

Component Alteration of a grtensor object: 

Applying routine expand to object Jump [K(dn, up) , f surf ] 
Applying routine simplify [trig] to object Jump [K(dn, up) ,f surf] 

CPU Time = .267 

> grdisplay (_) ; 

For the ssurf spacetime : 
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Jump from defaultMetric — Mint 



Jump [K{dn, up), f surf = ^'^^^ + cos( 



Jump [K{dn, up),fsurf]^ ' 



a( r ) sin( X ) 
Vcos(X)^ + cos(X; 



a( r ) sin( X ) 

We're nearly there, but Maple will not collapse e.g. ^/x^ to x unless it is sure 
X is real, or explicitly told to do so. We tell it to go ahead by using the routine 
simplify [sqrt , symbolic] . 

> gralter(_,sqrt) ; 

Component Alteration of a grtensor object: 

Applying routine simplify [sqrt] to object Jump [K(dn, up) ,f surf ] 

CPU Time = .066 

So if we take the FRW solution to be dust then we can establish that [gij] = 
and 



completing the junction of Schwarzschild to FRW. 



> grdisplay(_) ; 

For the ssurf spacetime : 

Jump from defaultMetric — Mint 

Jump[K{dn, up),fsurf] = All components are zero 
With the FRW solution restricted to uniform dust then we can establish that 
[gij] = and Kf = and hence a boundary surface exists provided we take 
R{t) = a^r) sinX and m = — R(t)'^ d'^ R{t) / dr'^ in Schwarzschild. 



4.2 Shells in Spherically Symmetric Static Space- 
times 

The evolution of thin shells in spherically symmetric spacetimes has been widely 
studied. Here we demonstrate how GRJunction is employed to determine the 
evolution of a thin-shell separating two spherically symmetric static spacetimes. 
This result contains as special case the classic analysis of Israel [^. 

We take M+ and M" as 

dsl = drl/f{r^) +rldQ^ - f{r_)dtl, (45) 
dsl = dr\/F{r+)+rldn^ - F{r+)dtl (46) 

where dVL^ = dO^ + s\n9'^d<jP' {9 and cj) are continuous through the surface). We 
define S in by r+ = R{t) and in M~ as r_ = R{t). On S we choose co- 
ordinates {9,(f),T). We seek to determine the equation governing the evolution 
of the surface and the stress-energy of the surface. 
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We demonstrate how to achieve this in the session below. Once again ev- 
erything from this point to the end of the subsection is Maple output. The 
specification of the surfaces follows exactly as in the Schwarzschild case in the 
previous example and so we omit the input portion of this process in the inter- 
ests of brevity. 

(Prior to this point we defined a surface sout in the metric ( |4^ ) which we 
labeled staticF and a surface sint in (^5|), labeled staticf .) 
We now identify these two surfaces by using join. 

> join(sout , sint) ; 

sout and sint are now joined. 
The default metric name is sout. 
The exterior metric is: sout 
The interior metric is: sint 

CPU Time = .050 
For the sout spacetime : 
Jump from defaultMetric — Mint 
Jump[g(dn, dn), sint] = All components are zero 
staticF 

We now calculate (the junction package object S3(dn,up)). Since this is 
non-zero it is clear that there is a thin shell seperating staticF and staticf. 

> grcalc(S3(dn,up)) ; 

Created definition for K(dn,up) 

CPU Time = 3.150 

> gralter(S3(dn,up) , factor); 
Component Alteration of a grtensor object: 

Applying routine factor to object S3(dn,up) 

CPU Time = .167 

To improve the appearance of the output we make use of GRTensor's ability 
to represent derivatives as subscripts (so dR{T)/dT — > R^) via the autoAlias 
command. This command resides in the grtools library, which we now load. 

> readlib (grtools) : 

Now we apply autoAlias via grmap. 

> grmap (S3 (dn, up) , autoAlias, 'x'); 
Applying routine autoAlias to S3(dn,up) 

> grdisplay(S3(dn,up) ) ; 

For the sout spacetime : 
Intrinsic stress — energy 
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Se' = -^{-2^iiR{T)) + Rr^F{RiT))-2,Jf{RiT)) + Rr^Rr^ 
+ 2^F(R(r)) + Vf(R(T)) + 2^F(R(r))+ 

- ^i{R{T)) + Rr^R{T)B{F){R{T)) - 2^f{R{T)) + Rr^Rir) Rr,r 

+ 7F(R(r)) + i?,2j^(r)D(/)(R(r)) + 2^F(R(r)) + VR(T)i?.,.) /(R(t) 
^F(R(r)) + V^f(R(r)) + V7r) 

5/ = -^(-2Vf(R(T)) + i?,2F(R(r))-2^f(R(T))+i?,2i?,2 
+ 2^F(R(r)) + i?^2f(j^(^-)) + 2^F(R(r))+i?^2^^2 

- ^Ji{R{r))+Rr^R{T)B{F){R{T))-2.Ji{R{T))+Rr^R{T)Rr,r 

+ ^F{R{t)) + Rr^R{T)B{f){R{T)) + 2^F{R{t)) + Rr^Rir) Rr,r) /(R(t) 
^F(R(r)) + i?^2 <^f(R(^)) + i?^2^^ 

Sr^ = -\{- ^JiiR{T))+Rr^F{RiT))-^f{RiT)) + Rr^Rr'' + ^FiR{T)) + Rr^f{RiT)) 

+ ^JF{R{T))+Rr^Rr'') / (r(t) ^Jf{R{t)) + Rr^ .Ji{R{T)) + Rr^ TT^ 

Note that since /(r) and -F(r) are unspecified their derivatives cannot be eval- 
uated. This results in the use of the Maple D derivative. We will demonstrate 
how to specify specific functions below. 

Next we consider the equation for the history of the shell HGeqn and the first 
integral of this equation evint 1 (which in the case of spherical symmetry stems 
from an identity). 

> grcalc(HGeqn, evlntl) ; 
Created definition for n(up) 
Created definition for R(dn,dn,up,up) 

CPU Time = 4.767 
Before displaying the results we factor the expressions. 

> gralter(_,f actor) ; 

Component Alteration of a grtensor object: 

Applying routine factor to object HGeqn 
Applying routine factor to object evlntl 

CPU Time = .317 

> grdisplay (_) ; 

For the sout spacetime : 
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History equation 



HGeqn = Q [- R{t) a{9,cP,r) ^f{R{T)) + Rr^B{F){Rir)) 

-2R{T)a{e,^,T)^f{R{T))+Rr^Rr,r-2R{T)a{e,^,T)^F{R{T))+Rr^Rr,r 

-R(r)a(0,</),r)^F(R(r)) + i?,2D(/)(R(r)) 

+ 4P{9,cl),T)^i{R{T)) + Rr^F{R{r)) + 4P{9,(P,r)yJi{R{T)) + Rr^Rr^ 

+ AF{e,<P,T) ^F{R{t)) + Rr^f{R{T)) + AFie,^,T) ^jF{RiT)) + Rr^ Rr^) /( 
R(r)v/F(R(r)) + i?.2 r( ^ ) ) + ^^^2^ ^ 

R(r)D(F)(R(r)) + F(R(T))-R(r)D(/)(R(T))-f(R(r)) \ 

R(t)2 ; 



e^;/n« = fi^^^ = — (- QAir^ a{9,(t>,Tf R(t)^ + r'^Firf - 2r'^F(r)f{r) + f(rf 
V 64 V 

+ 64(r2)3/2 7r2a(0,<^,rf R(r)3-32(r2)3/2 n"^ a{9,(l),Tf Rirf F{r) 
- 32 ( )3/2 a( 9, r )2 R( r )3 f( r ) + 256 tt^ (7( 0, (f), r )^ R( r )6 r^) /(tt^ 

a{9,cl>,rfRiT)\^fj 

Now we define functions for / and F. We consider the Israel thin-shell example 
and hence define /(r) = 1 and F{r) = 1 — 2m/r so we have a Minkowski interior 
and a Schwarzschild exterior. These functions are defined as Maple procedures. 
(This ensures that the defered derivatives will evaluate. Note that a simple 
substitution of e.g. f{R{t)) = 1 — 2m/R{t) would not affect D(f) (R(t)) since 
f (R(t) ) does not appear explicitly). A Maple procedure is declared in a state- 
ment of the form: procedureName := proc (arguments) procedure_body end: 

> f := proc(r) RETURN(l); end: 

> F := proc(r) RETURN ( l-2*m/r) ; end: 

To repeat Israel's dust shell analysis we must set the pressure to zero, which 
we now do. 

> grmap(_, subs, P(theta,phi ,tau)=0, 'xO ; 
Applying routine subs to HGeqn 
Applying routine subs to evlntl 



> gralter (_, power, expand) ; 

Component Alteration of a grtensor object: 

Applying routine simplify [power] to object HGeqn 
Applying routine simplify [power] to object evlntl 
Applying routine expand to object HGeqn 
Applying routine expand to object evlntl 
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CPU Time = 1.116 



> grdisplay(_) ; 

For the sout spacetime : 
History equation 

HGeqn = 
( 

1 a{6,(j),T)m 1 a{6,(j),T) Rt^t 1 a{9, (j),T) Rt^-, 



(1 771^ ( ^"^/^ 777, 

= -1 + T77 o ^ — + ,1^, X +47r^o-(6',0,r)^R(T' 
16 7r2 cr( 6*, r )2 R( T )*^ r3R(r) 

The above expressions correspond to the results given in Note that we did 
not have to integrate to obtain the first integral of the evolution equation. 

4.3 A complicated Minkowski junction 

In this section we describe the junction of the Minkowski metric in the form 

ds^ = '^dr^ + + + - dt^ (47) 

to a metric ds'^ of the same form, but with coordinates {R, U,^,T) and param- 
eter A. These metrics arise from setting m = in a Kerr metric and choosing 
u = a cos 9). The coordinate names r and R while standard are quite mis- 
leading since surfaces of constant "radius" describe spheroids with oblateness 
governed by a or A. We seek to join a surface of constant "radius" in M~ 
to M^. The surface in will be some function of R and U. This example 
demonstrates GRJunction's ability to handle non-spherical surfaces in a case 
where we know a boundary surface must result. 



The transformations from e.g. ( |47| ) to the Minkowski metric in spherical 
form (with coordinates f and 9) are given by 

u = a cos 9 (48) 
= r"^ + sin 6*^ 
fcos9 = rcos9. (49) 

If we take the definition of T," as r = X {X a constant) then we can use (^) 
to determine that the the definition of the corresponding surface in has 
parametric equation: 

= X - (R^ + A'^ -U^ - a"^ + ^(-2R'^A'^ + 2a'^A'^ -lA"^ + 2A'^U'^ + 
^ 4 
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2iR^A^ - 2R^A''a^ + 2R'A'' - 2R^A^U^ + a* A'' - 2a^A'^ + 
2a^A^U^ + ^8 - 2A^U'' + A^U^ + AR^U^a^A^)^'/^^)/A^)^'/^^ (50) 

and the relations = (,^*) are 

R = (-^2^ A.(-XV + a2^2^aV-a4 + (XV-2XV^2 

-2X2a"n2 + 2X2a6 + a'A"" + 2a^A^u^ - 2a^ A^ + a^u" - 2a^u^ + 
+AX^u^A^a'f'^) + X2 + a2 - ^,2)1/2 (5^) 

U = ^(-XV + 0^.42 ^aV-a4 + (xV-2XV^2_2^2^4^2 
av2 

+2X2a6 + a^A" + 2a'' AW - 2a^A^ + 

_2aV + a« + 4XV^2^2)i/2)i/2 ^53) 
$ = (53) 

t = r+(T). (54) 

In this case it not clear which sign we should choose for the normal vector in 
M+. However since we know a priori that a boundary surface must result if 
the package does not reach this result then we can consider the other choice of 
normal sign for n^. (In this case we choose the plus sign in (^). 

The package determines [gij\ = and that K~- is: 



(a2 - n2)VX2 + „2 



With careful direction of the computer simplification the junction package de- 
termines that is also given by (^5|). EJ 



4.4 Joining Kerr to Kerr 

In this section we consider joining the "interior" region of one Kerr spacetime 
in Boyer-Lindquist form 

ds^ = p (^^ + del^+{rl+a^) sine Jd<pl-dtl 

_^ 2mr_ (-^g-^g_2^^_ _ ^^_^2 
P 

p = r'^ + cos 5 = — 2mr_ + a'^ (56) 



^In an earlier attempt to determine Kf- (using Boyer-Lindquist coordinates) we encountered the 
Maple error "Object too Large" which occurs on 32-bit machines when an expression contains more 
than 64,000 terms. On present day workstations this limit can be reached in a matter of minutes 
and often little can be done to circumvent this Maple limitation. 
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with mass m and angular momentum a to the "exterior" region of another Kerr 
spacetime with mass M and angular momentum A (we use coordinates 
) ^± ) '/'it ) ^± ) for -M^). Such a problem would arise if a thin shell of matter 
was constructed around a Kerr black hole (i.e. a Dyson sphere). The general 
problem is an extremely difficult one as the "toy" problem m = M = indi- 
cates (see section |4.3| ) . In this section we limit the discussion to an expansion 
to first order in a and A, since to this order the surface in Ai^ is spherical. 

In the treatment of timelike spherical surfaces it is customary to use trans- 
formations such as r = R{t) and t = T{t) and then use Uau'^ = — 1 to eliminate 
dT{T)/dT. Note that in these spherical cases this produces grr = — 1 on S as 
desired. In all spherical cases T will be strictly a function of r but this fails to 
be true in Kerr. For Kerr spacetimes to achieve grr = — 1 we can use the same 
idea but if we blindly use t = T{t) we discover that in actuality t = T{t, 6). We 
can try again - using 6 as an argument and we do get tau as the proper time 
on the shell but now a 3T(r, 6)/ 89 appears in g^^, g^^ and g^^ . Since we know 

dT/dr is merely a function of 9 we can integrate trivially but this makes those 
components with a dT/d9 depend linearly on r and the metric components are 
explicitly dependent on proper time. Hence forcing r to be the proper time on 
the surface comes at considerable expense. Fortunately these problems do not 
arise in the order (a, A) expansions. 

To facilitate the matching of the gij on S we eliminate the g^^ terms on S 
by using a transformation to the zero angular momentum (ZAMO) frame in 
the definition of the surface. The transformations are 

r± = R,9± = 9, <P± = 4>-^± T±{t, 9),t±= T±{t, 9) (57) 

where VL = g(f,t/g(f>^- Using GRj unction we calculate the metric and second 
fundamental form on the surface and only then expand to order a, identify 
and calculate Sij. 

The package first determines 

ds"^ |s±= R^dP + sm9^d4>'^ - dr^ (58) 
and consequently [gij] = 0. 

The resulting stress-energy of the shell is (from this point on we drop the 
tildes on 9 and cp): 

f(R - M) - F{R - m) 
SirR'^fF 

R^sin^9 si 
3(am - AM) 



St 



c<P 

QT 

St 
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f{R-2M) - F{R-2m) 
AnR'^fF 



where / = y^l - 2m/R, F = y^l - 2M/R. 

Note the appearance of mixed term S'T which precludes a standard phe- 



nomenological interpretation (i.e. a as given by (23) is not an eigenvector of 



Sij). To first order we have 

^ am ^ AM 

and if we require r2_ = then the mixed terms vanish. This now permits 
standard phenomenological interpretation of the shell and we can now interpret 
the density and pressures in the usual way. Note that in this case the pressure 
is isotropic and hoop stresses do not arise. 

We can easily extend this calculation to the dynamic case with R R{t) 
in the above. The resulting stress energy tensor is: 



SI 



R^R{f -F) + f{R - M) - F{R - m) 



8-KR'^fF 
St = Si 
SI = R^sin'^e St 

3 AMf{R^{lfR + R- 2M))V2 _ amF{R^{R^R + R- 2m))i/2 



Sf 
SI 



Stt R'^fF 
f{R - 2M) - F{R - 2m) 



47ri?2/F 
where ' = d/dr. 



5 Verification 

In developing the junction package we have re-executed a number of calcula- 
tions performed in the literature to verify that the package can reproduce these 
results. We list some of the verifications we have performed in Table 5 (Maple 
worksheets for these verifications are available by ftp Q.) 

In addition we have verified the identities ( p^ and (|l8|) for a variety of 
spacetimes including a number of non-spherical wormhole spacetimes. 
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Problem 


Verified 


Comments 


Dust shell 0] 


evolution equation 


Minkowski-Shell-Schw. 
dynamics 


spherical 
symmetry Q 




Checked form of Kij for general 
spherical symmetry 


inhomogeneous 
slab cosmology 16] 


[9ij] = 
[Kij] = 


planar symmetry 


Schwarzschild 
wormhole [17| 


Sij, cr, P 


"Thin-shell" wormhole 

formed by joining Schw. exterior 

to Schw. exterior 


Schw- deSitter 
shell g 


Sij 


Spherical symmetry but uses 
coordinates {r,v,6,(j)) 


Rotating dust 
shell lis] 


[gij] = 0, K^^j to order 
Sij to order a 


match to Kerr exterior 
in small a limit 
(non-spherical surface) 


collapsing, rotating 
dust shell 


Sij 


Match to Kerr to order a 
(dynamic shell) 



Table 5: Results verified with GRj unction 
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